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We examine the conformations of a model for under- and 
overwound DNA. The molecule is represented as a cylindri- 
cally symmetric elastic string subjected to a stretching force 
and to constraints corresponding to a specification of the link 
number. We derive a fundamental relation between the Eu- 
ler angles that describe the curve and the topological linking 
number. Analytical expressions for the spatial configurations 
of the molecule in the infinite length limit are obtained. A 
unique configuration minimizes the energy for a given set of 
physical conditions. An elastic model incorporating thermal 
fluctuations provides excellent agreement with experimental 
results on the plectonemic transition. 
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I. HISTORY AND INTRODUCTION 

Conformations of a slender elastic rod were originally 
viewed as an interesting problem in classical elasticity 
theory. Kirchhoff jy was the first to make significant 
headway towards a complete solution. Almost a century 
later, as polymers became the subject of intense study, 
interest in the problem picked up once again Jl3| . After 
the discovery of biological polymers — nucleic acids and 
proteins — researchers recognized the importance of pre- 
dicting the elastic shape of linear molecules. The shape 
( "tertiary structure" ) of DNA and RNA plays an impor- 
tant role in the processes of replication and transcription. 
Because of this a number of authors have analyzed var- 
ious aspects of elastic DNA conformation [5,j$| for both 
closed (circular) ||f7|]8| and open (linear) [2| configura- 
tions. The approaches taken include Lagrangian mechan- 
ics J9|,|],@]|] i (numer ical) m olecular dynamics |@| and sta- 
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tistical mechanics plyL^ . Despite significant progress 
[p| , the equilibrium configurations of infinitely long open 
DNA have not been analytically described. Our main 
aims are to set up a formalism for obtaining equilib- 
rium configurations; to find one such conformation for 
stretched twisted DNA, and to set up a model of plec- 
tonemic transition to compare with experimental results 
& 

The work reported here is based on two major results 
obtained over a century apart. Kirchhoff jjj provided 
the basic framework of elasticity theory. He observed 
that the equations of equilibrium that describe an elas- 



tic rod are formally identical to the equations of mo- 
tion of a heavy symmetric top with one point fixed (see 
figure Q). The rod Hamiltonian is identical to the top 
Lagrangian, with arclength mapping onto time. This du- 
ality provides an added insight into the nature of the 
solutions. White's Theorem [Q provides another crucial 
analytical tool. This important theorem relates the link- 
ing number — a natural topological invariant of strings 
and ribbons — to two manageable components, each cal- 
culable in terms of locally defined quantities. According 
to White's theorem, Ck — Tw + Wr. Ck is the link- 
ing number, while Tw, the twist, monitors the twist of 
the molecule about its axis, and Wr, the writhe, records 
the contortions of the axis. This theorem greatly sim- 
plifies the problem of formulating a constraint on the 
linking number. In addition we make use of the results 
of Fixman and Kovac [^0| and Marko and Siggia jll 24 
to develop thermal effects of the plectonemic transition. 



II. THE MODEL 

The elastic model of DNA represents the molecule as a 
slender cylindrical elastic rod. To model external forces 
and torques the rod is stretched (in the z direction) by 
a force F and is required to have a fixed Ck (see fig- 
ure ^|top). The rod is parametrized by arclengh, s. At 
each point s we describe the rod by relating the local 
coordinate frame C to the frame Cq rigidly embedded in 
the curve in its relaxed configuration. The relationship 
between the stressed and unstressed local frames is spec- 
ified by Euler angles 9(s), 4>(s),ip(s) needed to rotate Co 
into C (see figure || bottom) . The shape of the backbone 
r(s) is traced out by the unit tangent t(s). n(s), a unit 
normal, keeps track of the twist, Tw. In this paper we 
will often omit s-dependence for brevity. We also make 
use of the notation x = 4-{x). Then, 

r(s) = f t(s') ds' (1) 
Jo 

t(s) — (sin 9 sin <fi, sin 9 cos <fi, cos 9) (2) 
n(s) = (cos 4> cos if) — cos 9 sin <f> sin tp, 

— cos (f) sin ip — cos 9 cos <j> sin ip, sin 9 sin ip) (3) 

Let the elastic constants of bending and torsional stiff- 
ness be denoted, respectively, by A and C, and let L be 
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the lengh of the rod. The energy of the twisted, stressed 
rod is the sum of bending and twisting energies and the 
potential energy produced by the stretching force F, Us- 
ing (|) and (§): 

Etot = E e i — F COS 9 = Etend + E tw i s t — F COS 9 = 

1 / (t) 2 ds + ^- ( ((n x h) ■ tf ds + 



(4) 



- Fcos6» 
Using (|J) and (^) in (0) we obtain 
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Etot = J ds '-j (j> 2 sin 2 9 + 2 ) + 

— (<j}cos9 + ^y -Fcos9. (5) 

The feature that sets this work apart from previous 
attempts and allows us to unambiguously determine a 
unique configuration for a given set of initial conditions 
is the constraint of maintaining a fixed linking number, 
Ck. Although linking number is usually associated with 
closed curves, the bound ends of our string allow us to 
define a fractional linking number for it. A caveat is 
that the local expressions we derive are only valid for the 
configuration (extended) we consider. White's theorem 
^ allows us to express Ck in terms of its components, 
Tw and Wr. Using © and(|): 



Tw — — 

2vr 



Wr = — 

477 



1 

(n x n) • t ds = 



i cos 9 



ds I ds 

!! JO 
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ipds (6) 
(7) 



A local expression for twist follows straightforwardly, as 
evidenced by the far right hand side of Eq. (||). The 
writhe, however, is not yet suitable for use as a Lagrange 
multiplier. To express it as an integral of a local quantity 
we use a theorem by Fuller || . The theorem allows us to 
define Wr locally using a diffcomorphism onto a reference 
curve, C r . (Once again, we must stress that a different 
configuration, e.g. a circular plasmid, would require a 
different reference curve, producing slightly modified lo- 
cal expressions). The writhe is expressed as 



Wr = Wr r 
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1 + t r • t 



ds , 
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where Wr r is the writhe of the reference curve. Finding 
a suitable reference curve proves crucial. The best choice 
is also the simplest - a straight line C r — (0,0, s). This 
gives t r = (0,0,1) and Wr r — 0. Substituting (g) into 
(p|) we obtain 



Wr = 



1 
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(cos# — 1) ds 



(9) 



Combining and (^) we are led to the simple expression 
for Ck: 



ds 



(10) 



Thus we have derived a simple conservation law that ex- 
presses the invariance of Ck. Inserting pdCk into the 
right hand side of Eq. (^|), with p a Lagrange multiplier, 
the expression to be minimized becomes 

H= [ — U 2 sin 2 9 + 9 2 
Jo 2 V 

+ — (j)cos9 + - F cosd - p(j) + ij^ ds (11) 

DNA conformations of minimum energy are found among 
the extrema of H. 



III. SOLUTIONS 

We find the extrema of TL by applying standard varia- 
tional techniques to (11). The resulting Euler - Lagrange 
equations for 6(s), <p(s) and ip(s) are: 



p(l — cos#) • p p (1 — cos 9) cos 9 
Asm 2 9 ^ = C Asln 2 ! ( ' 
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(13) 



A central goal is to find a unique conformation of the 
rod given a set of externally imposed constraints, F and 
Ck/L, where L is the length of the rod. We find that 
equations ( Jl3| ) and (|l2|) support two types of solutions. 
The first is a family of twisted vertical lines: 

9 = 0; = 0; ip(s) = {2irCk/L) s (14) 

The energy of the straight line follows directly from (||): 

2CT 



En 



(nCk/LY 



(15) 



The second family of solutions can be extracted from 
([TJ|) . Multiplying Eq. (|l^) by sin 2 9 and integrating once 
we can rewrite (O) in the following form (u = cos 9): 



ds = 



du 



^{l-u) + l{E Q -Fu){l-v?) 
du 



^ (u — a) (it — b) (u — c) 



(a < b < c) (16) 



These "writhing" solutions are characterized by the roots 
{a, 6, c} of the cubic polynomial in the denominator of 
( |l6| ) (see figure ||). One of the roots, either a or b, is 1. 
If u = 1 is a single root, then the configurations form 
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"superhelices." If u = 1 is a double root (a = b = 1) the 
molecule supports a soliton-like excitation (See figure ||). 
The quantity u = cos 6 takes on values between c and 
b, the quadrature turning points. The expressions for </>, 
ip and other quantities of interest follow by quadratures. 
The integrals are easily evaluated in terms of elliptic func- 
tions. 

We have investigated the properties of solutions to ( |l2] ) 
and (|l3|), with the bending and torsional stiffness, A and 
C, appropriate to DNA JT^,0] . Using both numerical and 
analytical methods we find that, given a particular Ck/T, 
the member of the writhing family with the lowest energy 
is the soliton configuration (a = b = 1) (see figures || and 
^) . The shape of this solution is defined by the following 
relationship between the arclengh s and u = cos 9: 



when the Ck/L approaches the critical value? Our results 
strongly indicate that the molecule attempts to loop over 
and pass through itself to shed a unit of Ck/L and thus 
starts to form a plectonemic bubble. In this sense the 
twisted rod is in a metastable state. The plectoneme 
plays the role of the "bounce" via which a system tun- 
nels out of the false vacuum. Beyond the transition to 
local instability the plectonemes ought to proliferate. To 
explore this scenario we formulate a very simple model 
of the plectonemic transition of stretched twisted DNA 
and compare its predictions with recent beautiful exper- 
iments by Strick et al. Jl9[ 



A. plectonemic transition model 



s(u) 



2F 



A(l 



■hi 



(17) 



where the lower root, c, is given by 
2 (nCCk/L) 2 



4AF 



AF 



1, 



(18) 



and 4> and tp are similarly determined. A very interesting 
quantity is the energy of the soliton and its relationship to 
that of the twisted line solution (which are both infinite 
in the limit L — > oo): 



E f 



olila 



Eu 



AE with AE 



AF r- 



— arctan 




> 



(19) 



It's clear from ( |l9| ) that the E so ii ton , while smaller than 
the energies of the other members of the "writhing" fam- 
ily, is always greater (by a finite amount) than the energy 
of the twisting line configuration satisfying the same con- 
ditions. 



IV. PLECTONEMIC TRANSITION 

Thus we find that an extended solution that minimizes 
energy and has a specified Ck/L it is always a twisted 
straight line. To check stability we perturb the straight 
line solution 6(s) = + 89(s). The perturbation calcula- 
tion shows that the non-trivial zero-energy mode satisfies 
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F 



o 



(20) 



Thus for Ck/L < j--k\/~F the straight line 6 — is sta- 
ble to small fluctuations. What happens to the molecule 



The model is diagrammed in fig. g[ In the follow- 
ing all quantities are normalized by the lengh L. DNA 
researchers prefer to use a = ACk to t/Cko to measure 
topological properties of DNA. Here Ck is the natural 
link of the unstressed DNA molecule; B-DNA has one 
right-handed twist every h = 3Anm. We will follow this 
notation. 

The molecule is constrained to have a total Link 
Cktot = &tot/h. The plectonemic fraction takes up X 
leaving 1 — X straight. The plectoneme has a radius R 
and a pitch P. The straight portion is twisted to its 
critical value dTw = h\J AF / (2nC) = ai/h. The actual 
twist is slightly below critical ]25|], but numerical results 
indicate that the precise value (which depends weakly 
on L) is adequately approximated by that of an infinite 
string [23]. Guided by "twist conservation" implied by 
eqns. (12) we assign the same rate of dTw to the plec- 
toneme. The remaining link, Cktot — {Tw p i + Cki), is 
absorbed by the plectoneme's Wr p i. Let us account the 
link distribution: 



Cki = Twk = -rX and Ck„ = Ck to t — Cki = 
h 



X = Wr p +Tw p 



(21) 



Because the plectoneme has the same rate of twist as the 
line, we can read off its writhe from eqn. (|2l|). At the 
same time the writhe of a plectoneme can be expressed as 
a function of P and R pl] , p2| . This gives us a constraint: 



Otot -oi ,.. XP 

= Yvr r , = 

h p 2n(R 2 + P 2 ) 



(22) 



Up to now we have not considered any thermal effects 
or corrections. Our aim is to build a formalism of obtain- 
ing equilibrium zero-temperature solutions about which a 
thermodynamic theory can be obtained (e.g. by consider- 
ing fluctuations) ijuj . However, because the experiments 
we are examining contain a regime in which thermal ef- 
fects play a significant role |19|we must consider them. 
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Marko and Siggia |2JJ,|ll|]have derived the free energy of 
a plectoneme in their examination of fluctuations about 
helical structures. To within order unity constants 



EL 



pl 



R 



(R 2 + P 2 ) 
((i?/r )- 12 + (vrF/ro)- 12 
A -i/3 ( i? -2/3 + (7rP) -2/3 



C (2iraiY 
2 V h 



(23) 



/ro 



The first two terms in ( |23| ) are elastic contributions from 
curvature and twist, respectively. The next line is the 
hard core interaction (ro ~ 1.75nm [M]). The last term 
is the entropic penalty incurred for winding too tightly 
p| . (It is interesting to note that although we include 
the last term in our model, its value is always negligible.) 
Setting the plectonemic fraction X we use (^2|) and ( |2^ ) 
to minimize E p i with respect to R and P. 

Next let us determine the thermal behavior of the 
straight-line segment. Such behavior for the untwisted 
rod has been examined in some detail by Fixman and 
Kovak pQj . Eqn. |2(] allows us to make use of their re- 
sults provided we replace F with F' = F — fr. Siggia 
provided a valuable summary of their results in an ap- 
proximate interpolation formula |24|,[ll| . We employ the 
above substitution in Siggia's result to solve for the ther- 
mal shortening of the straight portion of the molecule. Z 
is the actual observed extension: 



F 



4.4 
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= ((1 - Z/Xy 2 - 1 )+Z/X (24) 



k b T 

In the final analysis we compute the optimum plec- 
tonemic fraction X and the extension Z for a given Ck 
and F. The results are plotted in fig. || side by side with 
experimental results Jl9| . Because our model is a very 
simple one, and we have made no attemts to compute 
exact parameters (i.e. 'critical winding', etc.) we cannot 
claim perfect agreement. Nonetheless the resemblance is 
striking. Our model shows the shift from purely ther- 
mal behavior for very small a to a transition completely 
driven by elastic considerations for larger ct's and forces. 



FIG. 2. The elastic strand is described by a local coordi- 
nate frame, C. At each point C is related to an unstressed 
frame, Co, by Euler angles, The molecule is stretched by F 
and is required to have a constant Ck. 



FIG. 3. The behavior of the solutions is determined by the 
polynomial / (u) — 2F/A (u — a) (u — b) (u — c). The solid 
curve shows / (u) for a typical set of parameters. Motion is 
defined where f (u) > in the range u- < u < u+. The 
shaded curve shows a case where u + > 1. The dotted curve 
(u+ — 1) shows the set of parameters that minimize the en- 
ergy given a fixed Ck. This is the soliton configuration dis- 
cussed toward the end of the paper. 



FIG. 4. The writhing family of solutions. The solution with 
the lowest energy is the soliton. In the infinite-length limit 
the soliton and the twisted line have the same energy. 



FIG. 5. The extended and plectonemic phase coexist in the 
molecule. The plectonemic phase takes up a portion X. 



FIG. 6. A comparison of our predictions and the 
data of Strick et al. The families of curves are 
g = 0.102,0.043, 0.031, 0.001, andO from top to bottom. The 
stretching of the untwisted (a — 0) line is purely entropic; 
a = 0.102 transition is dominated by elastic energy. No at- 
tempts have been made to fit the data. 
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FIG. 1. A mapping of arclength onto time renders the 
equations governing symmetric elastic rods and spinning tops 
identical. The top left figure represents our elastic model of 
DN A (see figure ^) . The lower right figure is the conventional 
way of representing the motion of the top. The locus of the 
angles of inclination of the axis is represented as a curve on 
the unit sphere. 



4 



[1] KirchhoffG. J.F.Math (Crelle) 50, 285-313 (1859). 
[2] Benham C.J. Biopolymers 22, 2477-2495 (1983). 
[3] Fuller F.B. Proc. Natl. Acad. Sci. USA 75 No.8, 3357- 
3561 (1978). 

[4] White J.H. Am. J. Math. 91 No.3, 693-727 (1969). 
[5] Crick F.H.C. Proc. Natl. Acad. Set. USA 73 No.8, 2639- 
2643 (1976). 

[6] Schlick T. and Olsen W. Science 257, 1110-1115 (1992). 
[7] Tsuru H. and Wadati M. Biopolymers 25 2083-2096 
(1986). 

[8] TanakaF. and Takahashi H. J. Chera. Phys. 83(11) 6017- 
6026 (1985). 

[9] Le Bret M. Biopolymers 23, 1835-1867 (1984). 
[10] Gradshteyn, I. S. and Ryzhik, I. M. Table of Integrals, 
Series and Products, 5th Ed. (Academic Press, New 
York, 1980) 

[11] Marko J.F. and Siggia E.D. Phys. Rev. E 52, 2912-2938 
(1995). Marko and Siggia consider the thermal properties 
of a particular class of solutions - helices and plectonemes. 

[12] Hunt N. and Hearst J. J. Chera. Phys. 95(12), 9329-9336 
(1991). 

[13] Vologodskii A.V. and Cozzarelli N.R. Ann. Rev. of Bio- 
phys. and Biomol. Struct. V23:609-643 (1994). 

[14] vologodskii, A.V. et al, J. Mol. Biol. 227, 1224 (1992). 

[15] Shore D. and Baldwin, R.L J. Mol. Biol. 170, 957-981 
(1983). 

[16] Coleman, Sydney Aspects of Symmetry (Cambridge 

University Press, New York, 1985). 
[17] Boles T.C., White J.H. and Cozzarelli N.R. J. Mol. Biol. 

213, 931-951 (1990). 
[18] J. D. Gunton M. San Miguel and P. S. Sahni in Phase 

Transitions and Critical Phenomena, C. Domb and 

J. L. Lebowitz, eds. (Academic Press, New York, 1987), 

vol. 8. 

[19] Strick T.R., Allemand J.-F., Bensimon A., Croquette V. 
Science 271 1835-1837, (1996). 

[20] Fixman M. and Kovac J. J. Chera. Phys. 58, 1564 (1973). 

[21] Marco J.F. and Siggia E.D. Science 265 506-508, (1994). 

[22] White J.H. Mathematical Methods for DNA Se- 
quences, Ch. 9, CRC Press, 1989. 

[23] B. Fain, unpublished. 

[24] Bustamante C, Marko J.F., Siggia E.D., Smith S.B. Sci- 

ence265 1599 (1994). 
[25] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. A 38, 

364 (1988). 



5 



